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Two-dimensional magnetic garnets exhibit complex and fascinating magnetic domain structures, 
like stripes, labyrinths, cells and mixed states of stripes and cells. These patterns do change in 
a reversible way when the intensity of an externally applied magnetic field is varied. The main 
objective of this contribution is to present the results of a model that yields a rich pattern structure 
that closely resembles what is observed experimentally. Our model is a generalized two-dimensional 
Ising-like spin-one Hamiltonian with long-range interactions, which also incorporates anisotropy and 
Zeeman terms. The model is studied numerically, by means of Monte Carlo simulations. Changing 
the model parameters stripes, labyrinth and/or cellular domain structures are generated. For a 
variety of cases we display the patterns, determine the average size of the domains, the ordering 
transition temperature, specific heat, magnetic susceptibility and hysteresis cycle. Finally, we ex- 
amine the reversibility of the pattern evolution under variations of the applied magnetic field. The 
results we obtain are in good qualitative agreement with experiment. 

PACS numbers; 2001 PACS numbers; 75.70.Rf, 89.75.Kd and 61.43.Bn 
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INTRODUCTION 

Magnetic domains in magnetic garnet films are inter- 
esting in themselves because of the rich and complex 
pattern formation they display; stripes, cells, labyrinths, 
ordered hexagonal lattices, disordered froth, pentagonal 
structures along the front between the ordered lattice and 
the disordered froth, mixed states of stripes and cells, are 
some of the features observed experimentally [0, |[ |[ || . 
In addition, these patterns are quite universal, since they 
are similar to those observed in soap froth and in the an- 
nealing process of polycrystalline materials ||^, |^ . More- 
over, the transitions are driven by stress induced by an 
applied magnetic bias field H, a handy control parame- 
ter, which plays a role analogous to temperature and/or 
pressure in the melting transition of solids. 

Originally the structure of magnetic domains in ultra- 
thin magnetic garnet films attracted attention due to 
their potential application as magnetic bubble memo- 
ries [0. In fact, size and thermal stability of the domains 
are crucial to optimize recording performance. Due to the 
interesting physics related to these technological prob- 
lems, significant basic research efforts have been carried 
out lately [0, |[ The most extensively investigated 
systems are single crystal ferromagnetic thin films of the 
FeYBiGdGaO garnet type. They are 5 to 7 /im thick 
and exhibit a considerable anisotropy along the direc- 
tion perpendicular to the film, so that the magnetization 
remains along this axis except within the Bloch walls 
that separate opposite direction domains. These Bloch 
walls are of the order of 0.1 /xm thick. Different types 
of stable and metastable domain structures are gener- 
ated in such materials 0, ^, 0, depending on 



the transient magnetic fields applied during the thermo- 
magnetic cycling of the sample. Eye-appealing examples 
of the most common structures observed experimentally 
were obtained by Albuquerque and are found in refer- 
ence [g^ . The most prototypical examples, displayed 
in Figs. l|, can be grouped into several categories: par- 
allel stripes [Fig. |l|(a)], mixed states of labyrinths and 
bubbles [Fig. |l|(b)], pure labyrinths [Fig. 0(c)], and cell 
structures with a strong labyrinth component, as illus- 
trated in Fig. n 




FIG. 1: Stable domain structures observed experimentally by 
Albuquerque in a magnetic garnet. Obtained by varying 
the combination of perpendicular to the film and in plane 
applied magnetic fields. 

Often the bubble structure looks like soap froth, which 
originated the denomination "magnetic foam" . These 
magnetic foams are generated by cooling the sample, ini- 



2 



tially at a temperature T larger than the critical one Tc, 
to T < Tc in the presence of a weak magnetic field per- 
pendicular to the plane of the film. The magnetic bubble 
patterns result from a balance between the demagnetiz- 
ing energy (like long range dipolar interactions) and the 
energy required to create Bloch walls. This balance orig- 
inates an equilibrium domain size, which evolves as a 
function of the strength of the applied magnetic field. 
When this external field is increased, the bubble struc- 
ture coarsens and large bubbles, surrounded by smaller 
ones are observed, while the total number of bubbles de- 
creases. 

Another intriguing feature of these foams is the mem- 
ory effect. It consists in a reversible change of the foam 
topology under a cyclic variation of the applied magnetic 
field. For example, the bubble configuration observed 
under an applied magnetic field evolves into a labyrinth 
pattern at zero field, but the system returns to the pre- 
vious macroscopic configuration when the magnetic field 
is restored. It is the stability of the threefold vertices 
that seems to be responsible for this interesting memory 
effect |,|. 

Over the years significant effort has been devoted 
to the understanding of this rich and complex behav- 
ior [0, 0. For the particular case of magnetic froth Weaire 
et al. proposed a description based on an analogy 
with two-dimensional soap froth, including surface and 
compressibility energy terms. Sampaio et al. proposed 
a model Hamiltonian to describe the formation of mag- 
netic domains in strongly anisotropic magnetic thin films, 
on the basis of a ferromagnetic Ising Hamiltonian, includ- 
ing long-range (antiferromagnetic) dipolc-dipolc interac- 
tions and an external magnetic field. 

In this contribution we put forward an alternative 
model, addressing the problem of the magnetic domain 
pattern formation by means of a generalized spin Hamil- 
tonian. More precisely, we consider a three-state spin 
one Ising Hamiltonian, with competing long-range inter- 
actions, and we also include the anisotropy energy and 
an external field. 

This paper is organized as follows: after this Introduc- 
tion, in Sec. we present the model and the simulation 
technique. In Sec. the results reached by implementing 
our model are presented. The paper is closed in Sec. , 
where we provide a summary and draw conclusions. 



HAMILTONIAN AND SIMULATION 



Hamiltonian 



We model the system by a generalized spin one Ising- 
like Hamiltonian on a square lattice. The interaction 
between nearest neighboring spins is ferromagnetic and 
long-range competing dipolar interactions, between the 



n nearest neighbor sites (n > 1), are incorporated. We 
also include the magnetic field and the anisotropy. Ana- 
lytically our model Hamiltonian is 
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where the parameters Jij are the exchange coupling con- 
stants between the spin on the z-site and its first n near- 
est neighboring spins, located on sites j. Thus, Jy > 
(Jij < 0) imply ferromagnetic (antiferromagnetic) cou- 
pling. The anisotropy constant is defined such that A > 
favors spin orientation perpendicular to the plane of the 
sample. H^ and Hj are the intensity of the longitudinal 
and transverse external magnetic fields, respectively, and 
the i and j summations are performed over N{N — 1) lat- 
tice sites, since the primed summation implies that the 
self-interaction terms i — j are excluded. The transverse 
magnetic field and the anisotropy play a similar role, but 
while A may be considered as a microscopic character- 
istic of the system, and Ht are externally controlled 
parameters. 

If one chooses for the coupling constants Jij the dipolar 
interaction a square-angled labyrinth is generated, as in 
reference Instead, we adopt a long-range oscillatory 
potential of the RKKY-type. This way, and although the 
system is not metallic, use of an RKKY-like expression 
provides a set of oscillating coupling constants that de- 
pends on a single parameter k. The analytical expression 
we propose for is given by 
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where — |r 



"j\ is the distance between sites i and 
j. In the numerical calculations we consider interactions 
up to the seven-th nearest neighbor and for convenience 
we adopt as our unit of energy, here and throughout, the 
first neighbor exchange constant F, defined by 



_ cos(fca) _ 



(3) 



where a is the lattice parameter. This type of coupling 
guarantees that the interactions between nearest neigh- 
bors is always ferromagnetic and of the same strength, 
while the magnitude and sign of the coupling with the 
rest of the neighbors depends on k. 

Thus, our Hamiltonian differs from previous propos- 
als, like the one studied by Sampaio et al. ^ in at least 
three important aspects. We adopt: i) spin 1, instead of 
spin spin i, so that the domain walls may adopt finite 
width, since spin i implies abrupt transitions from one 
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domain to its neighboring one; ii) a long-range interac- 
tion that is not purely antiferromagnetic and decays as 



and iii) magnetic anisotropy is incorporated. As it C — - — --^ ■ (4) 



will become clear hereafter, the inclusion of these three {ksT) 

effects leads to significant qualitative and quantitative Finally, the spin susceptibility x at zero field is ob- 

differences in the results, and allows for a better descrip- tained from the magnetization fluctuations, and given 

tion of experimental observations. ]jy 



Simulation procedure 



(5) 



All of our calculations were performed on a iV x 
lattice, with N — 200, using a conventional Monte Carlo 
(MC) procedure to update the lattice spins at various 
temperature, T, and applied field, H, values. One Monte 
Carlo step (MCS) consists in the successive updating of 
N'^ spins chosen at random, following the Metropolis pro- 
cedure. Temperatures and energies are measured in units 
of r = 1, and k in units of the inverse lattice parameter 

Two different processes were investigated: quenching 
and relaxation, (a) Quenching: In order to quench the 
system the lattice temperature T is raised to infinity (all 
spin-states generated at random) and at the first Monte 
Carlo step T is reduced to a very low value (T — 0.005), 
after which the system is allowed to evolve for 1000 MCS. 
This way the system is prevented from exploring most of 
the available energy landscape and quite probably it is 
quenched into a metastable minimum energy state, dif- 
ferent from the lowest minimum, (b) Relaxation: On 
the other hand, thermal relaxation is achieved by slowly 
varying the temperature from a large initial value (T = 5) 
to a final temperature T = 0, in intervals of AT — 0.05. 
To be certain that thermal relaxation is indeed achieved 
we run 1000 MCS after every temperature step; there- 
fore, for a typical thermal relaxation run a total of ~ 10^ 
MCS are required, which is equivalent to 4 x 10^ individ- 
ual spin flips. 

During these processes, and as will be described in de- 
tail in Sec. , several thermodynamic properties and pa- 
rameters of the system are monitored: the energy E, 
the magnetization M, and the average domain size A 
are evaluated during the quenching process; along the 
relaxation process the specific heat C and the magnetic 
susceptibility x were are also monitored, as a function 
of T. In the latter case we perform an average over 20 
different replicas of the system in order to decrease the 
thermal fluctuations. The domain size A was evaluated 
by means of a Hoshen-Kopelman |l^ type algorithm. 
This modified algorithm computes the cluster size iden- 
tifying all clusters, calculating their surface and degen- 
eracy, and then evaluating the average cluster size. The 
specific heat C is calculated by means of the energy fluc- 
tuations and is given by 



RESULTS 
Model parameters 

The first part of our work was devoted to search for 
an adequate set of parameters, that generates magnetic 
domain patterns like the ones observed in magnetic foam 
experiments |l^, and displayed in Figs. |l|. The pro- 
cedure we employ is the following: we adopt a value of 
k, that in turn determines the values of the exchange pa- 
rameters Jij up to the seven-th nearest neighbor, without 
including anisotropy {A = 0) and for zero applied mag- 
netic field (H = 0). Next, the system is quenched from 
T — > cx) to a very low T value (T = 0.005). At the end 
of the run we keep or discard the choice of parameters by 
visual inspection. After selecting some interesting values 
of k we varied the anisotropy A. The k values we finally 
chose, following the above outlined scheme, yield the long 
range exchange functions illustrated in Fig. 0. 
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FIG. 2: Long range RKKY-like dipolar potentials used in the 
simulations. 

All the configurations are obtained after 1000 MCS, as 
explained in Sec. . This number of steps is far larger than 
what is required to reach steady state. To stress this fact 
we display in Fig. ^ the domain size A versus the first 
100 MCS is displayed. It is clearly noticed that after the 
first 30 MCS the system reaches steady state values. It is 
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also apparent that stripes and labyrinths are significantly 
smaller than the average cell dimension. Along with A, 
the total energy E and the magnetization M were also 
computed, but are not displayed graphically here as their 
time variation is similar to that of A. Thus, all in all it is 
evident from the behavior of these quantities, and from 
visual inspection of Fig. ^, that the system reaches steady 
state (metastable) equilibrium after relatively few MCS. 
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FIG. 3: Average domain size A versus number of Monte-Carlo 
steps (MCS). 

The domain structures generated, after quenching, give 
rise to three different patterns: stripes, labyrinths and 
cells. In all cases negative values of the anisotropy A < 
were adopted, to favor the presence of domain walls 
{S = sites). These sites emerge preferentially between 
opposite magnetization domains, since they lower the 
energy cost of growing a domain wall. Fig. [i] exhibits 
snapshots of the three types of domain structures ob- 
tained after quenching. They correspond to the following 
sets of parameters: i) fc = 4.2 and A — —0.5 (stripes); 
ii) k = 0.8 and A = -0.1 (cells); and finally iii) k = 1.2 
and A — —0.8 (labyrinths). 
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FIG. 4: Snapshots of the three types of domain structure 
obtained through simulations, after quenching. The stripe 
patterns (left) correspond to the parameters k = 4.2 and A — 
-0.5 . The cells (center) to fe = 0.8 and A = -0.1 . Finally, 
the labyrinths (right) to A: = 1.2 and A — —0.8 . 



As stated in Sec. these patterns emerge from the 
competition between the ferromagnetic first neighbor ex- 
change interaction F > 0, and the long-range oscillat- 
ing interactions given by Eq. ^ For the stripe structure 
(fc — 4.2, Fig. ^) the second neighbor exchange con- 
stant is antiferromagnetic and comparable in magnitude 
with the nearest neighbor one, while the third and fourth 
are ferromagnetic, leading to the formation of stripes of 
alternating spin orientation. This way, two frustrated 
couplings between nearest neighbors are compensated by 
four pairs of satisfied second neighbor interactions. The 
labyrinth case {k = 1.2), illustrated in Fig. is also a 
consequence of competing interactions between the ferro- 
magnetic first neighbor coupling and the rest, all of which 
are now negative and smaller in magnitude. In this case 
stripe domains, one lattice parameter wide, do not op- 
timize the energy thus yielding a domain width larger 
than in the previous case. Finally, when cells are gen- 
erated (see Fig. ^), both the first and second neighbor 
interactions are ferromagnetic, and all the rest are an- 
tiferromagnetic and of much smaller magnitude. Conse- 
quently, the system approaches full ferromagnetism and a 
quenched cell structure results from a frozen metastable 
spinodal decomposition. Later we will see that this cell 
structure, as well as the typical eddy structures of the 
labyrinth case, both originate in the quenching process. 

Before going on, it is important to stress the fact that 
each of the configurations in Fig. ^ originates solely from 
different Hamiltonian parameters, since both Jij (which 
is a function of fc) and A do vary. On the other hand, 
in the experiments each domain pattern is a metastable 
configuration driven by an external applied bias field H, 
as mentioned in the introduction. The effects of H on the 
domain patterns and on other physical properties will be 
addressed in Sees, and . 



Characterization 

Our main objective is to contribute toward an un- 
derstanding of the physical mechanisms responsible for 
the many qualitatively different configurations that two- 
dimensional magnetic systems do adopt. Thus, in order 
to properly characterize and systematize our conclusions, 
several tools and physical parameters are used to illus- 
trate the results we obtain. Clearly, direct pictures are 
the most straightforward form of visualization, and con- 
sequently we start providing representative snapshots of 
some of the typical patterns that the system adopts. 

To describe the equilibrium state we computed the 
specific heat C, the magnetic susceptibility x, and the 
magnetization M, as a function of temperature T. In 
addition, in Sec. and hysteretic and memory effects, 
respectively, will be presented and discussed. 

As already mentioned snapshots provide the most di- 
rect perception and they allow to observe directly the 



FIG. 5: Snapshots that illustrate the results of simulating 
the relaxation process for the labyrinth case, starting from a 
large temperature T^l: atr = 2a cell-like structure (left) 
develops. The ordering process starts around T = 0.8, as 
illustrated in (center). Finally, an ordered stripe-like pattern 
is reached at T = 0.1 (right). 



evolution of the system. The patterns displayed in Figs. ^ 
were obtained during runs at various T values, starting 
from a high initial temperature, fc^T > F, as illustrated 
in Fig. where a cell-like structure develops. As T is 
lowered the ordering process is triggered when the critical 
temperature T = Tc is approached (Fig. ^jb corresponds 
to T = 0.8 and a labyrinth pattern does emerge). Finally, 
the low temperature (Fig. ^ illustrates the T = 0.1 re- 
sult) the minimum energy configuration is attained. As 
expected, the low T stripe configuration reached in this 
way exhibits a higher degree of order than the quenched 
one, illustrated in Fig. In other words, the eddy pat- 
terns of the quenched labyrinth case are almost washed 
out by thermal relaxation, evidencing that the latter cor- 
respond to a metastable configuration. The experimental 
counterpart shows a very similar behavior, as reported by 
Mino et al. ||, |l|l for the (Hoi.aTbo.eBii.a) Feg O12 com- 
pound, when the system was relaxed from the quenched 
configuration under an oscillating magnetic field. This 
provides support to the fact that our model constitutes 
a reasonable choice. Moreover, the cell case also corre- 
sponds to a metastable pattern, obtained by freezing the 
spinodal decomposition of a quasi- ferromagnet. This is 
confirmed by relaxing the system, which leads to a per- 
fect fully magnetized (M — 1) ferromagnet. 

The MC thermal relaxation was carried out reducing 
the temperature from T = 5 to T = (we recall that 
energy units of F = 1, as defined in Eq. |^, are used 
throughout). From the E versus T plot of Fig. I it is ev- 
ident that a phase transition takes place at the inflection 
point. It is apparent as well that the energy sequence is 
Eiabyrinth > Egtripe > E^eiu as is verified by direct in- 
spection of Fig. ||. However, the most conclusive evidence 
for the phase transition is the specific heat C, which ex- 
hibits a characteristic peak (Fig. |^) that constitutes the 
signature of the transition. 

We recall that this curve -and all the others obtained 
during the relaxation processes- correspond to an aver- 
age performed over 20 independent runs. This allows to 
appreciate the benefits of such a procedure to attenuate 




FIG. 6: Energy versus T plots for the three types of patterns. 
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FIG. 7: Specific heat versus T plots for the patterns of Fig. 



the fluctuations inherent to computations on a relatively 
small system, of only 200x200 sites, by comparison with 
our preliminary results obtained with only one run. 
It is observed that the ordering temperatures are of the 
same order of magnitude as the nearest neighbor interac- 
tion F, which by definition is the same for all the samples. 
However, the critical temperature Tc is larger for the cell 
configuration, which is mainly ferromagnetic, and lower 
for the stripe and labyrinth patterns. 

For completeness we have also calculated the magneti- 
zation fluctuations, i.e. the spin susceptibility x at zero 
field, and the corresponding results are displayed in Fig. ^ 
for the three main structures we discuss. It is quite ap- 
parent that stripes and labyrinths, which order mainly 
antiferromagnetically, have a relatively low susceptibil- 
ity, while cells show a value of x more than an order of 
magnitude larger, consistent with their tendency to fer- 
romagnetism. 

Also, the thermal evolution of the magnetization, M 
versus T, portrayed in Fig. ^, is an evidence that the 
stripe and labyrinth structures are antiferromagnetic, as 
implied by the low amplitude fluctuations close to M = 0, 
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FIG. 8: Susceptibility x versus T for the patterns of Fig. |^. 



with (M) w 0, in sharp contrast with the cell structure 
results which exhibit a significant ferromagnetic magne- 
tization. All the above provides a picture that is self- 
consistent and agrees qualitatively with experiment. 
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FIG. 9: Magnetization M versus T for the patterns of Fig. |[ 



Hysteresis cycle 

In order to characterize the system as completely as 
possible, we also computed the magnetic hysteresis cy- 
cles by performing MC simulations at temperatures T < 
Tc 1, adopting the quenched state at zero applied field 
(H = 0) as the initial configuration, and thereafter in- 
creasing the magnetic field in increments of AH — 0.01 
(as always, in units of F). The magnitude of M was 
computed averaging over 1000 MCS for each value of the 
applied field intensity. This process was performed un- 
til saturation of the magnetization was achieved, after 
which the direction of the H field was reversed and the 
full hysteresis cycle computed. The results obtained are 
consistent with the previous discussion. The cell struc- 
ture exhibits a magnetic hysteresis typical of a soft ferro- 
magnetic system: that is, the system is easily saturated 



with a relatively low applied field (H « 0.05) and the 
remanent field is almost as large as the saturation field. 
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FIG. 10: Hysteresis loops (top) stripes for T — 0.2 and 0.5; 
(bottom) labyrinths for T = 0.05 . The virgin curves are 
marked with thick lines. 

In sharp contrast the stripe structure, illustrated in 
Fig. p^a, behaves quite differently. The hysteresis loop 
consists of two nearly rectangular portions, correspond- 
ing to positive and negative magnetizations M, shifted 
away relative to one another. The width of these rectan- 
gular portions shrinks considerably as the temperature is 
increased, but the relative displacement remains nearly 
constant. Thus, in spite of the fact that the equilibrium 
configuration is antiferromagnetic, at T = 0.2 the system 
remains saturated up to H r:! 1.3 implying that the en- 
ergy landscape is bumpy and with high energy barriers. 
In order to verify the preceding picture we also performed 
simulations at a higher temperature: T — 0.5. The cor- 



responding hysteresis cycle is also shown in Fig. 10a, and 
the dominant antiferromagnetic character is again evi- 
dent. The total magnetization has a sharp transition 
from zero to fully saturated, without going through inter- 
mediate values; as expected, the width of each of the rect- 
angular portions is significantly smaller than for T = 0.2 
and saturation is obtained for applied fields H ~ 0.8, i.e. 
not too different from the low temperature value. 

On the contrary the labyrinth configuration, plotted in 
Fig. where we represent the low temperature T — 
0.05 cycle, is rather peculiar. It requires large H values 
to reach saturation and the M versus H cycle is almost 
reversible {i.e. with very low remanence). Moreover, 
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two different regimes are clearly defined: for low fields 
the magnetization is essentially a linear function of H, 
and it is in this region where one can expect to obtain a 
memory effect. On the contrary, for large applied fields 
the system saturates, but when H is decreased it returns 
to the linear regime for the relatively large H « 2 field 
value. 



Memory effect 

Among the results obtained by Albuquerque per- 
haps the most striking one is the visualization of the 
memory effect illustrated in Figs. |ll|. ft consists in the re- 
versible change of the foam topology under a cyclic vari- 
ation of H. The "foam" observed under an applied field 
H in the top image of Fig. ^ evolves into a labyrinth 
pattern at zero field (middle figure), but the system re- 
covers the initial macroscopic configuration when H is 
restored to its original value, as seen in the bottom im- 
age of Fig. ft has been conjectured that it is the 
stability of the threefold vertices, that are so abundant 
in the upper and lower images displayed in Fig. 
is responsible for this interesting effect ||, ^ . 
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Unfortunately the memory effect is very difficult to 
simulate and our efforts in this direction have met with 
only modest success. In our quest we have investigated 
the labyrinth (fc = 1.2 and A ^ -0.8) and cell (fc = 0.8 
and A = —0.1) structures. In both cases we generated 
the initial structures by first quenching the system and 
afterwards applying a transverse magnetic field Ht = 0.5 
to favor the growth of a boundary (actually a "skin") 
of 5 = spins that separates the 5* = 1 from the S — 
— 1 stripes. It is important to stress that these Bloch 
walls (or "skins" ) assemble as a very narrow domain wall 
(one or two lattice parameters wide) but, in spite of their 
limited width, they are an important element to generate 
the memory effect. As a matter of fact, the effect fades 
away in the absence of these 5 = Bloch walls. However, 
it should also be remarked that, due to the length scale 
required to display our results in Figs. the "skins" are 
too narrow to be noticeable. 

Following the quenching and the application of H( a 
longitudinal field H^, lower than the saturation value, 
is also applied after which H is reduced to zero while 
performing 1000 or more MCS at every field value. The 
labyrinth structure thus obtained exhibits a sort of mem- 
ory effect, as can be seen in Figs, [l^a, ^2|b and 12;, ob- 
tained for Hi = 0, < 1.7 and = 0, respectively 
(Ht = 0.5 is kept constant during the cycle). After sev- 
eral cycles the structure does recover its overall original 
appearance, thus yielding a genuine memory effect. How- 
ever, there are differences with the experimental results: 
both the configuration we start from and the one we ar- 
rive at after cycling differ, both as far as the shape of 
the labyrinths and the width of the walls that separate 
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FIG. 11: Memory effect observed by Albuquerque ||l^ under 
cycling of the applied magnetic field. The top figure corre- 
sponds to H = 78.9 Oe, the center one to H = and the 
bottom figure again to H = 78.9 Oe. The close resemblance 
of the top and bottom figures is quite remarkable 



them is concerned, with experimental observation. On 
the other hand, in the cell structure case no memory ef- 
fect is observed, since the system is magnetically very 
soft. Thus, repeated cycling just yields a fully polarized 
lattice. 



SUMMARY AND CONCLUSION 

In this paper we model the magnetic relaxation and 
the formation of magnetic patterns in ultrathin films with 
anisotropy perpendicular to the film plane. To study this 
model the Monte-Carlo simulation technique was imple- 
mented in combination with a generalized two dimen- 
sional classical Ising-like Hamiltonian on a square lattice, 
including long range interactions, and transverse Hj and 
longitudinal Hi applied magnetic fields. For the dipolar 
interaction we adopted a long range oscillating coupling 
that extends all the way to the seven-th nearest neighbor. 

To present our results we start providing direct snap- 
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FIG. 12: Simulation results of the memory effect of Fig. |ll|for 
the labyrinth case (Ht = 0.5): (top) H; = 0, (center) EI; = 1.5, 
and (bottom) H; = . 

shots of the patterns that develop in the simulations. In 
addition, and to fully characterize the behavior of the 
system, we monitor the evolution of the total energy, 
the magnetization, the average domain size, the specific 
heat, the magnetic susceptibility and the hysteresis cy- 
cle. Altogether they provide a physical description which 
is qualitatively in agreement with experimental results, 
yielding stripe, labyrinth and cell patterns. The key pa- 
rameter that determines which type of pattern is gener- 
ated is the magnitude of k (the argument of the RKKY 
part of the potential given by Eq. ^ which modulates 
the oscillations and the relative weight of the long range 
interactions. On the other hand, the specific heat ver- 
sus temperature plot implies the existence of an order- 
ing temperature Tc below which pattern formation is ob- 



served after quenching the samples. 

Finally we address the striking memory effect observed 
by Albuquerque pO| , which consists in the reversible 
change of the foam topology under the cycling of the 
applied magnetic field. While our model fails to yield 
the shape and size of the domains observed experimen- 
tally, the patterns that evolve when the cycling process 
is implemented do display the experimentally established 
memory effect. 
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